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ABSTRACT 


Thomson  scattering  radiation  has  been  measured  in  terms  of  electron 
concentration  and  temperature  in  a  free  jet  facility  and  a  shock  tunnel  using 
a  1/2  joule,  Q-switched  ruby  laser.  A  detector  system  is  described  which 
attenuates  the  primary  ruby  wavelength  by  10^,  discriminates  between  ortho¬ 
gonal  polarization  states,  analyzes  the  radiation  spectrally  and  provides  a 
measure  of  radiation  background.  The  scattering  volume  was  approximately 
3  x  10-4  cm3  for  the  jet  facility  and  approximately  2  x  10"3  cm 3  for  the  shock 
tunnel.  A  temperature  range  of  3000  to  8000°K  and  an  electron  concentration 
range  5  x  10*4  to  5  x  10*5  cm~3  was  measured.  A  unique  method  of  measuring 
electron  temperature  is  described  using  a  passive  system  utilizing  a  retardation 
plate. 
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SECTION  I 
INTRODUCTION 


The  present  work  is  an  extension  of  the  laser  probe  program  initiated 
under  an  AEDC  contract  No.  A.  F.  40(600)-1147.  In  the  previous  work  a 
detector  system  was  built  which  permitted  the  determination  of  electron  den¬ 
sity  within  a  small  scattering  volume.  The  desirability  of  measuring  both 
electron  concentration  and  temperature  simultaneously  provided  the  motivation 
for  the  present  work. 

The  original  intent  of  the  program  was  to  measure,  during  a  laser  pulse, 
the  radiated  line  shape.  This  was  to  be  accomplished  by  electro-optical  tuning  - 
through  a  Fabry-Perot  interferometer;  however,  our  early  experiments  indicated 
that  this  approach  was  not  feasible  and  a  suitable  alternate  scheme  was  required. 
The  scheme  finally  used  which  best  fit  the  goals  of  the  contract  was  a  passive 
system  in  which  line  shape,  reflecting  electron  temperature,  could  be  obtained 
by  inference  from  the  ratio  of  two  separate  portions  of  the  line  shape.  In  addition, 
the  detector  system  had  to  retain  the  attenuation  characteristics  of  the  primary 
wavelength  of  the  original  detector  and  discriminate  between  the  radiation  back¬ 
ground  and  scattered  laser  radiation. 

The  detector  system  was  evaluated  in  both  an  arc  environment  and  a  shock 
tunnel  environment  in  which  simultaneous  measurements  of  electron.temperature 
and  concentration  were  made.  Although  the  arc  results  could  not  be  used  to 
validate  the  probe  due  to  high  impurity  levels,  the  shock  tunnel  measurements 
gave  good  agreement  with  equilibrium  theory.  Unfortunately  enough  tunnel 
measurements  could  not  be  made,  due  to  funding  limitation,  to  build  up  a  statis¬ 
tical  backlog  of  measurements  for  proper  evaluation. 
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SECTION  n 
DETECTOR  SYSTEM 


A.  General  Description 

The  detector  system  is  basically  a  differential  photometer  which  dis¬ 
criminates  the  incoming  radiation  in  terms  of  polarization  and  wavelength 
and  at  the  same  time  attenuates  the  primary  wavelength  to  about  5  orders  of 
magnitude.  The  optical  diagram  is  depicted  in  FIG  1  .  The  imaging  lens 
forms  a  1: 1  stigmatic  image  of  the  scattering  volume  on  the  entrance  slit  of 
the  detector.  The  slit  is  a  commercial  Hilger  bilateral  variable  slit  housing 
with  shutter  and  Hartmann  diaphragm.  The  collimating  lens  in  the  detector 
housing  collimates  the  radiation  onto  the  pair  of  spectral  mirrors  set  at  about 
7°  to  the  radiation.  The  spectral  mirrors  are  commercially  available  inter¬ 
ference  filters  without  blocking  filters  with  a  transmission  of  about  80%  and 
a  band  pass  of  about  5  A  centered  at  6943  A.  The  spectral  mirrors  are 
mounted  in  a  precision  goniometer  (see  Appendix  for  details)  which  can  be 
controlled  by  a  knob  from  outside  the  detector  housing.  The  primary  wave¬ 
length  (6943  A)  is  attenuated  by  the  spectral  mirrors  (10  traversals  between 
the  mirrors)  by  about  5  orders  of  magnitude  while  passing  radiation  5  A  away 
from  6943  A  with  a  high  optical  efficiency  (~75%).  The  field  lens  serves  to 
image  the  entrance  slit  onto  the  photomultipliers  and  to  collect  all  the  light 
coming  from  the  mirrors  (in  passing  through  the  spectral  mirrors  the  beam 
expands  to  about  1  inch  diameter).  A  broad  band  filter  ~100  A  and  80%  trans¬ 
mission  limits  the  radiation  passed  by  the  system.  The  Wallaston  prism  (1 
inch  diameter,  22  1/2°  divergence)  separates  the  incoming  radiation  into  two 
polarized  beams.  The  horizontal  beam  (parallel  to  the  plane  of  observation) 
is  used  to  measure  the  background  radiation  by  photomultiplier  #3  (PM  #3)  and 
does  not  contain  the  scattered  radiation  due  to  the  laser  light.  The  vertical 
polarization  containing  1/2  of  the  plasma  background  and  the  Thomson  scattered 
radiation  passes  through  a  100  ^  (6943)  retardation  plate  followed  by  a  second 
Wollaston  prism.  The  effect  of  the  retardation  plate  and  second  Wollaston  is 
to  wavelength  modulate  the  incoming  radiation  as' (sine) 2  and  (cosine)^  which 
falls  on  PM  #2  and  PM  #1  respectively.  All  the  lens  surfaces  were  anti¬ 
reflection  coated  to  give  an  overall  optical  efficiency  of  70%. 

In  operation,  the  detector  attenuates  the  primary  wavelength,  passes 
approximately  100  A  of  the  radiation,  presents  one-half  of  the  radiation  to 
PM  #3  not  containing  Thomson  scattered  light,  the  remaining  portion  being 
modulated  as  (cosine)2  and  (sine)^  and  presented  to  PM’s  #1  and  #2  respec¬ 
tively.  For  the  purpose  of  determining  electron  concentration  and  tempera - 
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ture,  the  sum  of  PM  #1  and  PM  #2  is  a  measure  of  the  concentration  and  the 
ratio  of  PM  #1  and  PM  #2  is  a  measure  of  the  electron  temperature.  PM  #3 
is  used  to  measure  the  magnitude  of  the  plasma  background  which  can  be  sub¬ 
tracted  from  PM's  #1  and  #2. 


To  facilitate  the  alignment  of  the  detector  system  with  the  scattering 
volume,  a  He/Ne  gas  laser  was  built  into  the  detector  assembly.  The  plug-in 
mirror  in  front  of  the  field  lens  is  used  to  direct  the  laser  beam  back  through 
the  spectral  mirrors,  collimator  and  entrance  slit  where  it  is  focused  in  space 
by  the  imaging  lens. 


R  Mirror  System  and  Retardation  Plate 

The  spectral  mirrors  (1/2”  x  2”)  are  regular  interference  filters  without 
the  blocking  filters.  These  are  mounted  in  a  precision  goniometer  which  is 
controlled  by  a  knob  outside  the  detector.  The  mirrors  were  designed  to  give 
minimum  reflectivity  (maximum  transmission)  for  6943  A  at  about  7°  angle  of 
incidence.  The  mirrors  were  tuned  initially  by  passing  light  from  a  tungsten 
lamp  through  the  mirror  system  and  observing  the  resultant  radiation  through 
a  2  meter  B&L  spectrogram.  The  angle  was  adjusted  until  there  was  maximum 
attenuation  for  6943  A.  In  addition,  the  mirror  holders  have  adjustments  for 
separation  between  mirrors  to  control  the  number  of  reflections  between  the 
mirrors.  In  the  final  configuration  the  separation  was  about  1"  giving  10  re¬ 
flections  betweep  the  mirror  pair.  The  spectral  sensitivity  of  the  mirror  pair 
due  to  finite  size  of  the  entrance  slit  is  about  262  A/radian.  Assuming  a  slit 
width  of  50  /a.  (used  in  these  experiments)  and  a  collimator  focal  length  of 
16  mm  amounts  to  ±1.  5  x  10"3  radians  or  ±0.  4  A.  This  corresponds  to  the 
approximate  line  width  of  the  burst  mode  of  the  laser. 


In  addition  to  tuning  the  spectral  mirrors,  the  retardation  plate  was  tuned 
in  a  similar  manner.  The  retardation  plate  is  a  quartz  birefugence  plate  ap¬ 
proximately  0.  75  cm  thick  and  1.  5  cm  diameter.  Using  a  tungsten  lamp  the 
retardation  plate  was  rotated  until  nulls  were  obtained  visually  through  the 
spectrograph.  Observations  were  made  on  both  the  (cosine)^  channel  and  the 
(sine)2  channel.  FIG.  2  shows  3  spectrograms  representing  the  resultant 
radiation  on  the  plane  of  the  PM's.  The  top  spectrogram  shows  the  resultant 
radiation  after  passing  through  the  spectral  mirrors.  The  second  shows  the 
wavelength  modulation  in  the  (cosine)2  channel  and  the  third  the  (sine)^  channel. 
The  lines  are  due  to  a  reference  neon  lamp.  These  spectrograms  were  taken 
without  the  broad  band  filter  in  place.  FIG  3  is  a  spectrogram  of  the  resultant 
radiation  through  the  (sine)2  channel  with  the  filter  in  place. 
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C.  Calibration  of  Detector  System 

In  order  to  successfully  use  the  detector  system  it  was  necessary  to 
provide  5  separate  calibrations  as  follows: 

(1)  attenuation  of  6943  A 

(2)  Transmission  curves  of  the  three  separate  channels 

(3)  Absolute  sensitivities  of  the  three  photomultipliers 

(4)  Overall  efficiency  of  the  three  channels 

(5)  Response  of  the  detector  system  to  differing  line 
shapes  corresponding  to  varying  electron  temperatures 

The  attenuation  measurements  were  performed  as  follows: 

The  laser ,operating  in  the  burst  mode, was  directed  to  a  glass  plate  at  an 
angle  of  incidence  of  about  15°.  The  light  is  reflected  into  the  detector  system 
with  the  detector  lens  focused  on  the  glass  plate.  The  energy  incident  on  the 
glass  plate  and  reflected  from  the  glass  plate  was  measured  with  the  Lear- 
Siegler  energy  monitor.  Approximately  50  millijoules  was  directed  into  the 
detector  system.  The  spectral  mirrors  were  adjusted  to  obtain  minimum  re¬ 
sponse  using  successive  firings  of  the  laser.  Due  to  the  inability  to  operate 
the  laser  in  the  Q-switched  mode  it  was  anticipated  that  the  attenuation  would 
not  be  near  theoretical  (-^lO-?)  due  to  the  excessive  line  width  (~.  3  to  .  4A). 
Measurements  were  made  varying  the  slit  width  of  the  detector  system  to  obtain 
a  measure  of  the  line  width.  At  20/4,  the  attenuation  is  less  than  10"®  and 
decreases  with  increasing  slit  widtn  to  about  10” 3.  To  test  the  hypothesis  that 
the  large  line  width  is  reducing  the  attenuation,  several  measurements  were 
made  with  the  laser  operating  near  threshold  (reduced  line  width).  It  is  found 
that  the  attenuation  coefficient  increased  by  a  factor  of  10.  In  the  final  testing 
of  the  detector  system  (in  an  air  arc)  using  a  Q-switched  laser  the  measured 
attenuation  was  less  than  10"5.  It  was  also  determined  that  due  to  scattering 


statistical  sampling  of  the  attenuation  coefficient  due  to  the  changing  line  shape 
of  each  firing  of  the  laser.- 

The  second  calibration  was  obtained  with  the  Heath  scanning  monochromator 
purchased  for  this  purpose.  The  detector  system  was  focused  on  the  exit  slit  of 
the  monochromator  which  was  fed  with  a  tungsten  ribbon  lamp.  The  transmission 
curves  (relative)  was  obtained  for  the  three  channels  representing  background 
light,  sin2  modulated  and  cos^  modulated  light.  Again  these  tests  were  conducted 
several  times  to  obtain  a  statistical  average  (FIGS  4  and  5  ). 
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The  absolute  sensitivities  were  obtained  in  the  following  manner.  The 
internal  mirror  which  normally  directs  the  internal  gas  laser  through  the 
spectral  mirrors  to  the  external  lens  of  the  detector  was  reversed  so  that  the 
laser  beam  is  directed  back  to  the  photomultipliers.  With  one  of  the  photo¬ 
multipliers  removed,  the  exit  beam  was  measured  with  a  calibrated  EG&G 
lite-mike.  A  calibrated  attenuator  for  6328A  (operating  wavelength  of  the 
laser)  was  then  placed  in  the  optical  path.  The  photomultiplier  was  then  re¬ 
placed  and  its  current  measured.  It  was  not  necessary  to  repeat  this  for  the 
other  photomultipliers  as  a  relative  calibration  was  made  by  replacing  .each 
tube  in  turn  in  the  same  channel  slot  with  the  detector  system  fed  with  a  tungsten 
ribbon  filament  lamp.  The  absolute  calibration  for  the  three  tubes  resulted  in 
the  following. 

PM  #1  (cos^  channel)  =  5. 40  x  10^  Q/cl  @  1200  volts 
PM  #2  (sin^  channel)  =  6.  66  x  10 15  0/cl  @  1200  volts 
PM  #3  (background  channel)  =  2. 29  x  10*6  gf/ cl  @  1200  volts 

where  0  is  photons  at  6943  and  cl  is  coulombs  of  charge.  In  addition,  a  relative 
calibration  of  the  three  tubes  was  made  at  1000,  1100,  1250  and  1300  volts  cathode, 
supply. 

In  order  to  measure  the  overall  efficiency  of  the  three  separate  channels, 
a  calibrated  blackbody  source  operating  at  900°C  was  used.  With  the  blackbody 
source  focused  on  the  slit  of  the  detector,  the  outputs  of  the  three  channels  were 
obtained  and  translated  into  photons/sec-A.  Using  the  relative  transmission 
curves  obtained  previously,  the  curves  were  integrated  (graphically)  to  obtain 
the  integrated  band  pass  for  each  channel  with  the  results: 

Channel  #1  =  59A 
Channel  #2  =  45A 
Channel  #3  =  80.  5A 

Using  these  results,  the  apparent  energy  to  the  detector  system  was  obtained 
and  then  ratioed  to  the  measured  outputs  with  the  following  results: 

Channel  #1  =48.  6% 

Channel  #2  =70.  9% 

Channel  #3  =  60.  7% 
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These  percentages  represent  the  peak  of  their  respective  transmission  curves. 

The  low  figure  of  Channel  #1  is  due  to  the  composite  of  the  rising  cos^  cycle 
and  the  declining  slope  due  to  the  spectral  mirrors.  Experimentally,  the  sum 
of  channels  #1  and  #2  equals  #3  indicating  very  little  overall  loss  of  light.  The 
overall  efficiency  (sum  of  all  three  channels)  was  measured  at  61%.  As  the 
interference  filter  at  the  peak  transmission  of  the  system  is  about  80%,  this 
indicates  an  overall  efficiency  of  75%.  As  the  retardation  plate  and  one  lens 
is  not  coated  in  the  system  combined  with  the  10  reflections  from  the  spectral 
mirrors,  this  figure  is  close  to  what  is  expected  (~80%). 

In  order  to  determine  the  response  of  channels  #1  and  #2  to  changing  line 
shape  (electron  temperature)  it  was  necessary  to  compute  the  theoretical  profiles 
for  equilibrium  air  conditions  and  integrate  the  theoretical  profiles  over  the  ex¬ 
perimentally  determined  transmission  curves  to  obtain  the  ratio  response  in 
terms  of  electron  temperature.  The  method  used  for  computation  of  the  theoretical 
profiles  was  described  in  our  previous  report  and  will  not  be  repeated  here.  For 
the  air  arc  experiments  the  resultant  calibration  curve  for  electron  temperature 
is  shown  in  FIG  6  .  As  indicated  before,  the  electron  concentration  is  determined 
by  the  sum  of  Channels  #1  and  #2  after  subtraction  of  background.  The  method 
used  will  be  described  under  the  section  on  air  arc  experiments. 

For  the  shock  tunnel  experiments  where  equilibrium  cannot  be  assumed, 
it  was  necessary  to  compute  the  ratio  of  Channels  #1  and  #2  as  function  of  Te, 
o*  and  Ne  where;  oi  is  the  l/KD  with  K  =  4jL/>»n  9/2.  and  D  the  Debye  length 
(see  previous  report).  * 


D.  Electronics  and  Read  Out 

The  signal  data  handling  system  design  was  based  on  the  desire  to  provide 
maximum  flexibility  and  to  fully  utilize  the  received  scattered  laser  light  signal. 
A  block  diagram  of  the  data  handling  system  is  shown  in  FIG  7  .  The  three 
photomultipliers  are  RCA  Type  4526  tubes,  selected  for  quantum  efficiency  of 
at  least  7%  at  the  ruby  laser  wavelength.  PM’s  #1  and  #2  receive  the  Thomson 
scattered  laser  light  as  well  as  any  plasma  background  light  which  is  passed  by 
the  detector  optical  system.  PM  #3  sees  only  background  light  and  is  used  to 
provide  a  measurement  of  this  background  which  can  be  used  to  subtract  out  the 
background  component  from  the  other  PM  outputs. 

The  gate  generator  is  a  tunnel  diode  junction  transistor  level  detector 
which  provides  a  constant  amplitude  gate  output  for  a  time  about  fifty  mano- 
seconds  longer  than  the  duration  of  the  laser  pulse  in  order  to  allow  for  the 
electron  transit  time  dealy  in  the  PM’s.  The  input  to  the  gate  generator  is 
provided  by  the  Laser  Energy  Monitor  whose  silicon  photodiode  generates  a 
signal  each  time  the  laser  pulses.  Emitter  followers  couple  the  gate  to  each 
of  the  gated  integrators. 
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PM  #1  and  PM  #2  outputs  are  each  returned  to  ground  through  a  fast 
recovery  hot-carrier  diode  (Hewlett-Packard  5082-2800  series)  and  a  fifty  ohm 
resistor  which  shunt  the  integrating  capacitor  at  the  input  to  an  FET  amplifier. 
When  the  laser  fires,  a  gate  pulse  is  developed  across  the  fifty  ohm  resistor 
to  reverse  bias  the  diode  and  divert  the  PM  output  into  the  integrating  capacitor. 
This  signal  is  amplified  by  about  30  db  and  coupled  into  an  FET  storage  amplifier  • 
for  presentation  on  the  display  oscilloscope.  The  storage  amplifier  is  of  a  sum¬ 
ming  type  such  that  the  PM  signals  obtained  during  each  laser  pulse  are  totalized. 
The  storage  time  constant  is  about  ten  milliseconds  in  order  to  minimize  the 
decay  between  pulses.  This  circuit  was  necessitated  by  the  nature  of  the  laser, 
which  multiple-pulsed  even  when  Q-switched.  The  output  signal  for  a  typical 
multiple  laser  pulse  is  shown  in  FIG  8  which  also  shows  the  laser  energy 
monitor  output  and  resultant  gate  for  a  single  pulse. 

Unfortunately,  while  this  data  handling  system  is  sound  in  principle  and 
behaved  well  when  tested  using  a  steady  background  signal,  it  proved  inadequate 
when  subjected  to  the  rapid  fluctuations  of  intensity  of  the  arc  plasma  in  the  facility 
used  to  evaluate  the  laser  probe.  These  fluctuations  prevented  the  establishment 
of  a  quiescent  voltage  level  at  the  amplifier  input  and  caused  sufficient  zero  shift 
in  the  storage  amplifier  output  to  make  meaningful  measurements  virtually  im¬ 
possible.  While  a  solution  to  these  problems  is  undoubtedly  attainable  by  redesign 
of  the  integrator  and  use  of  a  sample  and  hold  circuit  in  place  of  the  simple  storage 
circuit,  this  was  unnecessary  for  the  purpose  of  this  program.  Problems  with 
the  laser  led  to  the  use  of  a  different  one  which  did  not  multiply  Q-switch,  thus 
eliminating ; the  need  for  the  summing  storage  circuit.  Use  of  shorter  coupling 
time  constants  in  the  amplifier  minimized  the  effects  of  the  arc  fluctuations  and 
allowed  for  direct  display  of  the  amplified  PM  output  signal. 
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FIGURE  1.  OPTICAL  SCHEMATIC  OF  DETECTOR  SYSTEM 
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FIGURE  2.  SPECTROGRAMS  TAKEN  THROUGH  DETECTOR  SYSTEM  (SEE  TEXT) 
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FIGURE  4.  SPECTRAL  TRANSMISSION  CURVE  FOR  (SINE)2  CHANNEL 
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FIGURE  5.  SPECTRAL  TRANSMISSION  CURVE  FOR  (COSINE)2  CHANNEL 
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FIGURE  6.  CALIBRATION  CURVE  FOR  ELECTRON  TEMPERATURE 
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TOP:  GATE  GENERATOR  OUTPUT 

BOTTOM:  LASER  INTENSITY  MONITOR  OUTPUT 


STORAGE  AMPLIFIER  OUTPUT 

TOP:  0.  275  JOULES  OUTPUT  LASER 
BOTTOM:  0. 145  JOULES  OUTPUT  LASER 

FIGURE  8.  LASER  AND  ELECTRONICS  CIRCUIT  CHARACTERISTICS 
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THOMSON  SCATTERING  EXPERIMENT 


A.  Air  Arc  Experiments 

The  experimental  arrangement  was  similar  to  the  previous  work  except 
the  arc  generator  was  a  free  burning  jet  and  the  detector  axis  was  oriented  150° 
to  the  laser  beam.  The  electric  vector  of  the  laser  beam  (plane  of  polarization) 
was  oriented  perpendicular  to  the  plane  of  observation. 

I  _ 

1L  Gerdien  Free  Jet  Generator 

The  Gerdien  Free  Jet  arc  facility  (FIG  9)  is  one  of  a  family  of  plasma 
arcs  available  at  GE-RESO  which  is  based  upon  the  constructed  arc  principle. 
This  principal  is  shown  schematically  in  FIG  10.  The  sketch  shown  is  that 
of  a  tandem  unit.  It  is  composed  of  five  chambers,  two  electrode  chambers, 
two  vortex  chambers,  and  a  plenum.  Gas  enters  the  vortex  chambers  tan¬ 
gentially,  some  of  which  passes  into  the  rear  electrode  chamber,  the  re¬ 
maining  portion  passing  into  the  plenum.  That  portion  which  enters  the 
electrode  chamber  entrains  whatever  material  may  evolve  from  the  electrode 
and  the  gas  carries  this  with  it  as  it  is  exhausted  overboard.  The  gas  entering 
the  plenum  is  heated  as  it  passes  through  the  constrictor  and  again  in  the 
plenum  where  it  mixes  with  the  portion  of  heated  gas  entering  from  the  op¬ 
posite  side.  In  a  tunnel  operation  the  plenum  gas  is  expanded  to  lower  pres¬ 
sure  through  a  nozzle  while  in  free  jet  operation  the  gas  expands  only  sub- 
sonically. 

The  Gerdien  Free  Jet  operates  basically  on  the  same  principle  with  some 
modification  (FIG  11).  Only  one  side  of  this  unit  is  as  shown  in  FIG  10;  the 
other  electrode  is  cylindrical  and  the  arc  column  exiting  from  the  constrictor 
actually  arcs  to  the  cylindrical  wall.  By  using  an  external  field  the  arc  attach¬ 
ment  point  to  the  metal  surface  is  continually  moved  about  the  periphery 
thereby  maintaining  its  integrity.  The  gas  then  passes  through  this  second 
electrode  and  exits  from  a  subsonic  nozzle. 

i  Power  is  supplied  to  the  facility  from  a  500  KW  DC  mercury  rectifier  , 
ufoit.  Electrical  ballast  is  used  in  series  with  the  arc  to  vary  arc  power.  * 
This  is  provided  by  variable  length  of  water  cooled  stainless  steel  pipe. 
Although  ballast  can  be  changed  over  a  range  of  ohmage  values,  it  cannot 
be  varied  while  the  arc  is  in  operation. 
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For  the  present  study  test  flows  were  established  over  a  range  of 
arc  operating  conditions.  Maximum  arc  power  was  250  KW,  with  maximum 
current  being  1100  amps.  While  a  number  of  different  nozzles  geometries 
were  used,  all  the  data  were  generated  in  subsonic  flows;  observations 
were  made  of  the  exit  flow  from  an  orifice  which  was  1/8”  to  3/8"  in  diameter 
depending  upon  the  particular  test  These  spectrographic  measurements 
were  obtained  at  a  d' stance  of  12  mm  above  the  nozzle  exit 


•* 


2.  Laser  System 

The  laser  system  used  for  these  experiments  was  a  Lear-Siegler 
LS-140-Q  ruby  system  consisting  of  a  laser  head,  energy  storage  bank 
and  control  unit 

t 

The  laser  head  is  an  air  cooled  (fan)  dual- elliptical  cavity,  pumped  by 
linear  flash  lamps.  The  ruby  rod  lies  along  the  common  focus  of  the  two 
ellipses.  The  Q-switch  module  consists  of  uranyl  glass  separately  pumped 
by  a  linear  'flash  lamp.  One  end  of  the  optical  circuit  is  terminated  by  a 
resonant  reflector  with  the  output  end  of  the  ruby  rod  flat  and  unelevated. 

The  system  employs  three  energy  storage  units  rated  at  200  MFD  at 
4000  volts  delivering  1600  joules  each.  Consequently,  the  system  can  be 
charged  to, 4800  joules.  The  system  can  be  manually  or  at  1  pulse/minute. 

For  both  the  arc  and  tunnel  experiments,  the  laser  was  operated  at 
about  0. 7  joules  on  the  Q-switched  mode  with  a  pulse  width  of  30  nanoseconds. 

A  Lear-Siegler  laser  energy  monitor  (Ml-2)  was  used  in  conjunction  with 
the  laser.  This  is  a  sampling  device  in  which  beam  splitters  are  used  to 
sample  a  portion  of  the  laser  light  as  it  passes  through  the  monitor.  Light 
from  the  beam  splitter  illuminates  a  diffuse  spherical  reflector  which  is 
viewed  by  a  high  speed  silicon  photodiode  through  a  set  of  optical  attenuators. 

A  holding  feature  of  the  electrometer  readout  current  permits  readings  to  be 
made  directly  in  joules  within  thirty  seconds  after  firing.  In  addition,  the 
instantaneous  power  can  be  read  on  an  oscilloscope.  Both  the  power  and 
energy  measurements  are  considered  to  be  accurate  to  ±10%. 

In  the  Q-switched  mode,  the  beam  was  elliptical  in  shape  with  a  divergence 
of  2. 5  by  4  milliradians.  The  spectral  line  width  was  less  than  0. 1A.  In 
addition  to  the  main  laser  beam,  the  laser  head  produced  a  highly  diverging 
beam  of  scattered  light  which  required  some  baffling  in  front  of  the  focusing 
lens  in  spite  of  the  high  attenuation  in  the  detector  system.  A  glass  lens, 

2"  diameter  a  16"  focal  length  was  used  to  focus  the  beam  into  the  arc  jet. 

The  lens  formed  a  spot  approximately  2x3  mm. 


i 
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3.  Alignment  Procedure 

As  mentioned  above,  a  gas  He/Ne  laser  was  incorporated  into  the  detector 
system  to  facilitate  alignment  of  the  detector  with  the  laser.  A  piece  of  car¬ 
bon  paper  was  placed  on  the  jet  exit  orifice  and  the  laser  fired  leaving  a  hole 
in  the  carbon  paper.  The  ruby  laser  was  adjusted  until  the  beam  was  at  the 
desired  height  — approximately  12  mm  above  the  exit  orifice.  The  gas  laser 
in  the  detector  system  was  turned  on  and  the  plug-in  mirror  inserted. 

The  focused  beam  formed  by  the  gas  laser  was  used  to  adjust  the  detector 
to  the  mark  or  hole  left  on  the  carbon  paper  by  the  ruby  laser  beam.  In 
this  manner  it  was  possible  to  align  the  detector  to  within  0. 1  mm  of  the 
ruby  beam.  After  alignment,  the  plug-in  mirror  was  removed  and  the 
detector  closed.  Periodically,  the  alignment  was  checked. 

4.  Experimental  Procedure 

Due  to  the  fact  that  both  electron  temperature  and  concentration  measure¬ 
ments  are  made  during  a  single  firing  of  the  laser,  it  was  not  necessary  to 
calibrate  the  arc  provided  the  assumption  of  equilibrium  and  clean  flow  can 
be  made.  The  philosophy  of  this  work  was  that  no  matter  how  erratic  the 
arc  behaved,  the  simultaneous  measurement  of  temperature  and  concentration 
should  reflect  equilibrium  values  and  hence  validate  the  probe  measurements. 

With  the  air  flow  off,  spurious  scattering  was  found  to  be  orders  of  mag¬ 
nitude  greater  than  the  expected  signal.  This  was  attributed  to  Mie  scattering 
in  apparently  dirty  air.  Consequently  the  level  of  spurious  scattering  could 
not  be  obtained  during  arc  runs  to  subtract  from  the  Thomson  signal.  The 
assumption  was  made  the  spurious  scattering  did  not  contribute  appreciably 
to  the  Thomson  signal.  Occasionally,  a  large  signal  would  be  obtained  in¬ 
dicating  a  particle  or  particle  cloud  was  struck  and  the  results  disregarded. 

Periodically  a  check  was  made  on  the  experimental  data  for  distinguising 
features  of  Thomson  scattering: 

1.  Synchronization  with  the  laser  pulse 

2.  Polarization  in  the  same  sense  as  the  input  polarization 

3.  Misalignment  of  the  detector  with  the  laser  beam  to 
insure  that  the  recorded  signal  is  not  some  anaxnoly  of 
the  plasma. 
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After  the  arc  was  initiated,  the  laser  was  fired  after  the  arc  stabilized 
(**3  sec).  The  outputs  of  PM’s  #1,  #2  and  #3  were  recorded  on  scopes 
(see  electronics  section).  FIGS  12, 13,14 show  some  typical  oscillograms 
obtained  during  an  arc  run  from  PM’s  #1  and  #2.  All  runs  were  made  with 
the  PM  voltage  at  -1200  volts,  the  entrance  slit  of  the  detector  at  50 ^  and 
a  slit  height  of  1  mm.  The  average  energy  of  the  laser  was  recorded  for 
each  run  with  the  energy  monitor,  varying  between  0.  4  to  0. 7  joules.  The 
initial  delay  evidenced  in  the  oscillogram  is  due  to  the  transit  time  through 
the  optics  and  the  PM’s  as  the  trace  was  initiated  from  the  energy  monitor. 


5.  Data  Reduction  and  Results 

In  order  to  relate  the  observed  voltage  signals  from  PM’s  #1  and  #2 
into  electron  concentrrtion,  it  was  found  more  expedient  (rather  than  working 
from  a  normalized  calibration  curve)  to  derive  a  working  equation.  This  is 
derived  in  the  appendix  with  the  results, 

Ne  =  $T  (measured)  3>  60  x  1010  (cm*3)  (1) 

J  (measured) 

where  is  the  sum  of  photons  received  by  PM’s  1  and  2  and  J  is  the  energy 
measured  in  joules.  The  electron  temperature  was  obtained  by  ratioing  PM’s 
1  and  2  after  converting  to  photons  using  the  calibration  curve  FIG  6  . 

The  initial  experiments  utilizing  the  free  jet  showed  the  existence  of 
large  high  frequency  fluctuations  in  the  intensity  level  of  the  plasma  radiation. 
This  resulted  in  unacceptably  large  fluctuations  in  the  voltage  at  the  input  to 
the  gated  integrator.  These  were  amplified  and  coupled  to  the  output  causing 
a  base-line  shift  which  made  accurate  measurement  of  the  gated  integrator 
output  during  the  laser  pulse  difficult.  The  amplifier  circuits  were  therefore 
modified  to  have  a  20  KHZ  low  frequency  cut-off.  This  reduced  the  effect  of 
the  arc  fluctuations  to  an  acceptable  level.  FIG  15  is  a  typical  oscillogram 
of  the  arc  background  fluctuations  obtained  from  PM  #3  (background  channel). 
The  sweep  is  50/t sec/cm  and  amplitude  is  2  V/DIV.  The  output  of  PM  #3 
was  coupled  into  a  1  M-TL  resistor  paralleled  by  a  120  PF  capacitor. 

Approximately  150  air  arc  tests  were  made  in  an  effort  to  validate  the 
instrument  response.  Initially,  using  copper  and  carbon  electrodes  for  the 
arc,  measurements  of  electron  concentration  and  temperature  yielded  in¬ 
consistent  results.  The  ratio  response  which  is  a  measure  of  electron 
temperature  did  not  correlate  with  the  total  response  which  is  a  measure  of 
electron  concentration.  The  results  were  also  not  reproducible  yielding  a 
range  of  temperatures  and  concentrations.  The  results  also  indicated  that 
there  is  still  some  interference  from  solid  particles  in  the  arc  jet. 
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The  arc  configuration  was  changed  to  both  carbon  electrodes  which  pro¬ 
duces  a  higher  temperature  jet.  With  the  laser  probe  aligned  approximately 
12  mm  above  the  exit  orifice  of  the  jet,  runs  were  made  in  which  only  the 
electrical  ballast  to  the  arc  was  changed.  As  the  ballast  was  lowered,  the 
instrument  registered  higher  temperatures  as  expected;  however,  at  the 
lowest  ballast  conditions  an  unexpected  decrease  in  temperature  was  ob¬ 
served  (Table  1).  Due  to  the  uncertainties  experienced  with  the  arc  jet,  it 
was  decided  to  perform  a  spectrographic  profile  of  the  jet  at  the  position  of 
the'laser  scattering  experiments. 

A  2  meter  B&L  spectrograph  in  conjunction  with  a  dove  prism  was 
used  to  characterize  the  arc  jet.  The  dove  prism  was  necessary  to  invert 
the  jet  image  such  that  the  spectrographic  slit  cut  radially  across  the  jet 
at  12  mm  from  the  exit  orifice.  The  initial  spectrogram  showed  atomic 
copper  lines,  calcium  atomic  and  ion  lines  (impurities)  and  atomic  lines 
of  nitrogen  and  oxygen  in  addition  to  bands  systems  of  CN  and  C2.  Using 
a  calibrated  tungsten  strip  an  absolute  calibration  was  made.  The  copper 
lines  reflected  a  temperature  of  5000°K.  By  using  the  ratio  of  the  atomic 
to  ionic  lines  of  calcium,  which  is  sensitive  to  temperature  and  electron 
concentration,  an  electron  concentration  of  10^  cm “3  was  measured. 

The  equilibrium  value  is  6.  4  x  10^  indicating  that  impurities  may  be 
increasing  the  concentration  over  equilibrated  air.  This  deviation  was 
not  enough  to  account  for  the  large  deviation  measured  with  the  electron 
probe  (a  factor  of  10  over  equilibrium  conditions).  In  addition,  there  is 
some  doubt  as  to  the  radiating  regions  of  both  calcium  and  copper  (i.  e.  if 
calcium  is  radiating  from  a  different  volume  than  the  copper  then  the  above 
results  would  be  inconclusive).  Consequently,  an  absolute  line  intensity 
measurement  was  made  on  the  nitrogen  atomic  line  at  7424A.  In  addition 
using  a  standard  inversion  technique^),  it  was  possible  to  obtain  the  radial 
temperature  profile  of  the  arc  jet.  This  is  a  lengthy  procedure  and  subject 
to  many  errors;  however,  the  results  obtained  can  sometimes  be  useful. 

It  was  found  that  from  the  edge  of  the  nitrogen  radiating  region  (approxi¬ 
mately  7  mm  from  the  center  of  the  arc  jet)  to  the  center,  the  temperature 
ranged  from  7100°K  to  8000°K.  This  is  rather  a  flat  profile  and  quite  un¬ 
expected  for  a  jet  of  this  type  as  a  higher  temperature  gradient  would  be 
expected.  This  would  seem  to  suggest  that  the  jet  column  is  swinging  from 
side  to  side  in  some  manner  which  would  tend  to  reduce  the  observed  tem¬ 
perature  gradient.  If  this  be  the  case,  it  would  also  explain  the  large 
variation  of  results  obtained  with  the  laser  probe  as  during  a  typical  firing 
of  the  laser,  the  scattering  takes  place  at  different  locations  in  the  arc  jet. 


21 


AEOC-TR-71  -69 


To  resolve  the  large  variation  of  electron  concentration  and  temperature 
obtained,  high  speed  movies  of  the  arc  jet  were  taken  with  a  Dynafax  camera. 
This  camera  takes  pictures  of  1  second  duration  and  about  40 ja-  seconds 
apart.  Movies  were  taken  at  about  the  same  time  that  the  laser  is  fired  after 
the  arc  is  initiated  (~  3  sec).  The  movies  showed  sharp  local  variations  of 
luminosity  from  frame  to  frame  (FIG  16).  In  some  cases  a  given  point  in 
the  jet  column  shows  no  luminosity  while  at  later  times  shows  maximum 
luminosity  (complete  blackening  of  the  film).  As  the  luminosity  observed  is 
probably  due  to  the  CN  molecule  emission  an  attempt  was  made  to  obtain  high 
speed  movies  of  the  jet  using  an  oxygen  filter  (7947A)  with  the  camera.  This 
proved  unsuccessful  as  the  luminosity  was  not  sufficient  to  expose  the  infrared 
film,  even  at  reduced  exposure  times. 


As  the  laser  probe  results  were  obtained  with  a  sampling  time  of 
about  50  nanoseconds  while  the  dynafax  pictures  were  of  1 u.  second 
duration,  it  would  be  difficult  to  obtain  any  correlation  between  the 
two.  In  addition,  as  the  luminosity  of  the  CN  molecule  is  not  necessarily 
linearly  related  to  electron  temperature  or  concentration  we  were  not 
able  to  reach  any  conclusion  as  to  the  validity  of  the  laser  probe  results. 
Further,  to  obtain  enough  of  a  statistical  sampling  to  correlate  the  results 
with  the  spectrographic  results  would  require  thousands  of  runs  which  is 
not  feasible.  Consequently,  at  this  stage  of  the  program,  it  appears  that 
validation  of  the  laser  probe  is  not  possible  in  the  arc  jet  source  being 
used.  Validation  would  require  a  more  stable  and  uniform  source  than  is 
currently  available. 

Although  we  were  confident  that  the  probe  results  are  accurate  (based 
on  our  previous  calibrations),  our  inability  to  validate  the  results  dictated 
that  we  discontinue  our  testing  in  the  present  source.  Consequently,  the 
laser  probe  was  set  up  in  our  High  Enthalpy  shock  tunnel  which  is  described 
in  part  B. 

6.  Discussion  of  Results  of  Table  I 

Although  Fastax  movies  of  the  arc  jet  indicated  that  at  the  time  of  the 
laser  firing,  the  position  of  the  jet  with  respect  to  the  laser  and  detector  can 
vary  giving  a  range  of  temperature  and  concentrations,  it  was  noted  at 
about  6000°K,  the  electron  concentration  indicated  an  air  equilibrium  value. 
At  lower  temperatures  it  would  appear  that  impurities  control  the  electron 
temperature  which  is  reflected  in  the  very  high  electron  concentrations.  At 
the  lowest  ballast  (highest  current  density)  it  also  appears  that  so  much 
material  is  drawn  off  from  the  electrodes  that  the  plasma  temperature  is 
again  controlled  by  the  impurity  level.  Between  6000  and  7000°K  the  results 
seem  to  indicate  that  both  air  electrons  and  impurity  electrons  contribute  to 
the  temperature  which  is  reflected  in  the  higher  than  equilibrium  values  for 
concentration. 
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B.  Shock  Tunnel  Experiments 

1.  High  Enthalpy  Shock  Tunnel  Facility 

This  facility  has  a  high  performance  electrical  discharge  heated-helium 
driver.  The  driver  is  4  1/2  feet  in  length  and  normally  is  operated  with  an 
electrically  insulated  sleeve  having  a  2  1/2  inches  I.  D.  The  sleeve  has  an 
inner  coating  of  bonded  RTV  rubber,  which  provides  the  electrical  insulation 
required.  Operation  is  off  of  a  768  KJ  capacitor  bank  with  a  maximum  voltage 
ratio  of  40  KV.  The  driver  tube,  made  of  stainless  steel,  is  17  feet  long 
and  has  a  2  inches  L  D.  In  the  current  study,  the  facility  was  operated  as 
a  shock  tunnel  in  the  reflected  mode.  A  converging-diverging  conical  nozzle 
was  located  at  the  downstream  end  of  the  tube.  This  nozzle  has  a  l°50f 
half-angle.  Measurements  were  taken  at  the  3.  6  area  ratio  station  where 
the  nozzle  diameter  was  0.  38".  A  sketch  of  the  facility  is  given  in  FIG  17. 
Here  an  MHD  accelerator  is  located  at  the  downstream  end  of  the  driven 
tube  instead  of  the  conical  nozzle.  A  dump  tank  is  downstream  of  the  nozzle 
or  accelerator  section.  In  these  experiments,  the  MHD  accelerator  was 
not  used. 

2.  Experimental  Arrangement 

The  experimental  arrangement  used  in  the  high  enthalpy  shock  tunnel 
is  depicted  in  FIG  18.  The  detector  was  orientated  173°  to  the  laser  beam 
with  the  plane  of  polarization  perpendicular  to  the  plane  of  observation.  Due 
to  the  excessive  scattering  of  the  laser  beam  from  the  entrance  window,  it 
was  necessary  to  provide  a  light  dump  composed  of  a  carbon  block  with  a 
tapered  hole.  The  surface  reflection  from  the  focusing  lens  provided  enough 
laser  light  for  an  EG&G  lite-mike  which  was  used  to  indicate  the  time  of  the 
laser  pulse.  It  was  also  found  necessary  to  provide  a  1-meter  long  collimator 
tube  composed  of  1  cm  diameter  aperture  every  1”  along  the  collimator  tube. 
This  was  dictated  by  the  widely  divergence  scattered  beam  originating  within 
the  laser  cavity.  The  Lear-Siegler  energy  monitor  was  used  to  provide  the 
energy  measurements  of  the  laser  beam. 

Due  to  the  angle  at  which  scattering  was  observed,  the  scattering  volume 
was  an  ellipse  cylinder  1  mm  x  12  mm  x  50  u, .  The  experimental  results 
were  obtained  at  a  point  approximately  7  mm  downstream  of  the  flow  from 
the  exit  nozzle.  FIGS  2 6  and  27  are  photographs  of  the  experimental  set-up. 
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3.  Experimental  Procedure 

Initially,  the  detector  was  aligned  to  the  laser  beam  using  the  procedure 
outlined  in  the  arc  tests.  The  tunnel  wac;  pumped  down  and  the  laser  was 
fired  to  indicate  the  level  of  spurious  scattering.  Only  the  (cosine)2  channel 
indicated  any  spurious  scattering.  The  level  of  spurious  scattering  was 
constant  to  the  normalized  laser  output  indicating  only  apparatus  scattering 
rather  than  Mie  scattering.  During  the  pre-test  runs  the  Lite-Mike  was 
calibrated  in  terms  of  its  output  against  the  energy  monitor.  Typical  os¬ 
cillograms  taken  during  a  pre-test  are  shown  in  FIGS  21,22,23.  FIG  21 
shows  the  output  of  PM  #1  at  20^  sec/ CM  and  5  V/CM.  FIG  22  shows  the 
output  of  PM  #2  at  50yc<  sec/ CM  and  5  V/CM  and  FIG  23  shows  the  output 
of  PM  #3  (top  trace)  at  1  V/CM  and  5Qx  sec/ CM  and  the  output  of  the  Lite- 
Mike  (bottom  trace)  at  100  MV/CM  and  5Q^<  sec/ CM. 

During  the  tests  conducted  on  the  high  enthalpy  shock  tunnel  flow,  a 
pressure  gage  in  the  driven  tube  wall  was  used  to  provide  a  master  trigger 
signal  from  which  the  triggering  of  a1!  scopes  and  of  the  laser  flash  lamps 
was  obtained  through  appropriate  delays.  Two  oscilloscopes  were  used  to 
display  the  outputs  of  sidewall  pressure  gages  further  downstream  whose 
signals  were  displayed  in  pairs  on  each  scope  to  provide  a  measure  of 
shock  wave  velocity.  A  third  scope  was  used  to  display  the  reflected  region 
pressure  as  a  check  on  the  state  of  the  test  flow  at  the  time  of  laser  Q- 
switching.  Oscilloscope  £5  displayed  the  output  of  PM  #3  and  of  the  EG&G. 
Lite-Mike  energy  monitor  as  a  check  on  background  light  and  energy  in  each 
laser  pulse.  Scopes  #6  and  #7  displayed  the  amplified  outputs  of  PM's  #2 
and  #1,  respectively. 

In  general,  the  Q-switching  of  the  laser  occurred  about  900  microseconds 
after  triggering  of  the  flash  lamps.  Therefore,  a  delay  of  about  100  micro¬ 
seconds  after  the  main  trigger  would  assure  laser  Q-switching  during  estab¬ 
lished  test  gas  flow  in  the  nozzle. 
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4.  Data  Reduction  and  Results 

In  order  to  translate  the  PM  outputs  into  electron  temperature  and 
concentration,  it  was  necessary  to  reprogram  the  line  shape  equation  for 
&  =  173°.  As  constant  pressure  cannot  be  assumed  in  this  case  as  it  was 
in  the  free  jet,  the  program  consisted  of  solving  the  line  shape  equation 
for  constant  Os  and  a  range  of  temperatures  from  4000°K  to  8000°K  and 
computing  the  response  of  each  channel.  As  a  given  oi.  +  Te  with  0  con¬ 
stant  define  an  Ne,  it  was  possible  to  plot  Ne,  Te,  and  Os  against  p  the 
ratio  of  the  (sine)*  channel  to  the  (cosine)  2  channel.  This  is  reproduced 
in  FIG  24.  To  show  that  the  sum  of  channels  #1  and  #2  are  independent 
of  <*>  within  the  temperature  range,  the  results  of  the  computer  program 
provided  a  plot  of  Ne,  cK  ,  response  of  #1,  #2  and  the  sum  of  #1  and  #2 
against  0/joule.  This  is  reproduced  in  FIG  25.  The  program  listing  and 
the  results  of  the  program  is  reproduced  in  the  Appendix. 

In  order  to  facilitate  the  determination  of  the  electron  concentration  a 
working  equation  similar  to  the  arc  runs  was  derived  from  the  geometry  of 
the  arrangement.  The  calculations  are  relegated  to  the  Appendix  with  the 
result, 


Ne 


6. 80  x  109 


0  measured  sum 
Joule  (measured) 


With  the  determination  of  Ne  and  measured  ratio  it  was  then  possible  to 
determine  the  electron  temperature  from  FIG  24  independent  of  knowing 
the  pressure  or  state  of  equilibrium.  For  each  run  the  value  of  the  constant 
apparatus  scattering  was  subtracted  from  the  observed  output  before  the 
computation  was  made.  The  Ne  computation  can  be  made  from  calibration 
curves  in  FIG  25  or  from  the  above  equation. 

Due  to  limitation  of  funds  only  two  tunnel  runs  could  be  made.  The 
oscillograms  from  one  of  them,  Test  304  is  shown  in  FIGS  26,  27  and 
28.  In  FIG  26  the  upper  trace  shows  the  varying  response  from  PM  #3 
and  the  bottom  trace  the  output  from  the  energy  monitor.  As  can  be  seen, 
the  laser  Q-switched  twice.  The  first  pulse  produced  enough  scattering  to 
show  up  in  PM  #3.  FIG  27  indicates  that  PM  #1  saturated  ('■v  10  volts), 
consequently,  a  determination  of  Te  and  Ne  could  not  be  made.  However, 
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it  is  instructive  to  compute  at  what  level  PM  #3  would  show  a  laser  signal. 
With  a  resistive  load  of  100  Kil.  the  change  in  current  during  the  laser 
pulse  was  4  x  10-®  amps.  With  a  sensitivity  of  2.  3  x  10*6  0/ci  this  cor¬ 
responds  to  9.  2  x  10lu  0/sec.  Assuming  a  laser  pulse  duration  of  30  x 
10_»  secs. ,  PM  #3  received  2800  0.  The  energy  of  the  pulse  was  0.  44 
joules  and  hence  the  output  from  PM  #3  represents  =  6400  0/joule.  With 
a  transmission  of  0.  61  for  this  clannel  this  corresponds  to  ~10, 550  0/joule 
at  the  entrance  to  the  detector  within  the  band-pass.  The  output  from  PM  #2 
(FIG  35  )  indicates  an  energy  input  of  ^  108, 000  0/ joule  using  a  similar 
analysis.  Due  to  the  fact  that  PM  #1  contains  some  constant  spurious  scat¬ 
tering  plus  Thomson  signal,  the  output  of  PM  #1  was  calculated  for  the  con¬ 
ditions  of  the  run  at  the  time  of  the  laser  pulse  (i.  e.  Te  =  6600°K  and  Ne 
5  x  10 15).  Calculating  the  expected  return  for  PM  #1  and  adding  the  amount 
of  spurious  scattering  as  calculated  from  the  energy  monitor,  the  input  to 
PM  #1  at  the  entrance  to  the  detector  amounted  to  300, 000  0/joule.  As 
PM  #3  records  the  sum  of  PM’s  #1  and  #2,  this  indicates  that  PM  #3  sees 
approximately  2.  5%  of  the  Thomson  signal.  This  is  not  unexpected  as  the 
Wollaston  prism  separating  the  two  polarization  states  is  not  anti -reflection 
coated  and  Fresnel  reflection  from  the  prism  finds  its  way  into  PM  #3.  This 
is  not  particularly  objectionable  as  PM's  #1  and  #2  were  calibrated  in  the 
presence  of  the  surface  reflection. 

As  the  laser  Q-switched  twice,  it  was  possible  to  obtain  a  result  on  the 
second  pulse  (with  less  energy)  as  neither  PM  #1  or  PM  #2  were  saturated. 
The  second  pulse  in  PM  #2  is  not  visible  in  the  reproduction.  Inserting  the 
appropriate  values  the  total  return  to  PMrs  #1  and  #2  was  44, 000  0.  From 
the  working  equation  we  obtain  Ne  =  2. 8  x  1015.  The  ratio  of  #2/#l  was 
0.  6uu  and  hence  from  FIG  31  we  obtain  Te  =  6000°K.  An  error  statement 
is  discussed  in  the  Appendix. 

In  FIGS  29,  30,  31  the  outputs  are  reproduced  for  Test  #305.  In 
this  case  the  laser  Q-switched  three  times.  The  first  pulse  again  saturated 
PM  #1.  The  second  pulse  indicated  particle  scattering  as  the  ratio  obtained 
was  about  10:1  (no  PM  #2  output).  Only  pulse  number  3  could  be  used  for 
analyses.  For  this  pulse  PM  #1  recorded  43, 500  0  and  PM  #2  31, 900  0. 

The  ratio  0. 732  and  the  Ne  obtained  from  the  working  equation  indicated  an 
Te  =  6500OK  and  Ne  =  5. 7  x  10*5  were  indicated.  The  comparison  of  these 
results  with  the  value  measured  from  pressure  and  velocity  measurements 
are  indicated  below: 
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Test  #304 


Test  #305 


Ne  (Laser)-.  >»6000°K 
Ne  (Actual)  J(^6000oK 
Te  (Laser)  A*"*  2. 8  x  10l5(cm”3) 
Te  (Actual)/^  3  x  10 15  (cm “3) 


6500 °K 
6530°K  ie 
5. 7  x  10l5(cm-3) 
5.  5  x  10l®(cm“3) 


Although  two  results  are  not  an  indication  of  the  validity  of  the  probe,  the 
value  of  these  probe  results  are  described  in  a  paper  of  C.  J.  Harris  of  the 
Environmental  Sciences  Laboratory  staff  which  is  reproduced  in  the  Appendix. 
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FIGURE  9.  GERDIEN  FREE  JET  FACILITY 
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FIGURE  10.  SCHEMATIC  TANDEM  GERDIEN  ARC  HEATER 
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FIGURE  11.  SCHEMATIC  OF  GERDIEN  FREE  JET 
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FIGURE  12 
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FIGURE  15.  OSCILLOGRAM  OF  ARC  BACKGROUND 
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FIGURE  16.  PORTION  OF  FASTAX  MOVIES  OF  ARC  JET 
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TABLE  I 


Height  Above  Jet  -  12  mm 
Slit  Width  -  50/# 

Slit  Height  -  1  mm 


Ballast  (Ohms) 

0.  52 

0.  52 

0.  39 

0.39 

0.26 

0.26 

0. 195 

0. 195 

0. 195 

Energy  (Joules) 

0.  49 

0.54 

0.  54 

0.  52 

0.48 

0.46 

0.48 

0.52 

0.  51 

Run  No. 

117 

118 

119 

120 

121 

122 

123 

124 

125 

Vl  (MV) 

505 

540 

775 

425 

375 

350 

400 

475 

425 

V2  (MV) 

65 

95 

95 

100 

115 

125 

50 

50 

60 

01  (14.  5  0/MV) 

7,  300 

7,830 

11,240 

6, 160 

5, 440 

5,  080 

5,  800 

6,890 

6,160 

02  (11.8  0/MV) 

770 

1, 120 

1, 120 

1, 180 

1,  360 

1,480 

590 

590 

700 

(01  +  02) 

8,070 

8,950 

12,360 

7,340 

6, 800 

6,560 

6,  390 

7,480 

6,860 

R 

0. 105 

0.143 

0.  100 

0. 192 

0.250 

0.291 

0. 102 

0.  086 

0. 114 

Te  (OR) 

3,000 

3, 800 

2,950 

5,400 

6, 700 

7,  300 

2,950 

3000 

3, 200 

Nex  10-14(CM-3) 

5.9 

6.0 

8.2 

5.1 

5.1 

5.  1 

4.8 

5.2 

4.8 

Ne  x  10“14(CM-3)+eq. 

6.  4  x  10’4 

.04 

6  x  10"4 

0.9 

6.1 

8.8 

6  x  10-4 

6  x  10-4 

7  x  10 

*  Indicates  the  electron  concentration  at  the  measured  temperature  for  equilibrated  air, 
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FIGURE  18.  EXPERIMENTAL  ARRANGEMENT  FOR  PROBE 
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FIGURE  20.  PHOTOGRAPH  OF  LASER  PROBE  EXPERIMENTAL  SET-UP 
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FIGURE  21 
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FIGURE  24.  CALIBRATION  CURVE  FOR  TUNNEL  ELECTRON  TEMPERATURE 
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FIGURE  25.  CALIBRATION  CURVE  FOR  TUNNEL  N 
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FIGURE  26 
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FIGURE  29 
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FIGURE  30 
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FIGURE  31 
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APPENDIX  I 

CALCULATION  OF  WORKING  EQUATIONS 


In  order  to  facilitate  the  determination  of  the  electron  concentration  from 
the  raw  data,  it  is  convenient  to  write  the  expression  for  Ne  as, 


where  J 0  =  number  of  photons 

-^2-  =  solid  angle  of  detector 

6  eff.  =  effective  cross  section  for  Thomson  scattering 
Vscat.  =  scattering  volume  as  seen  by  the  detector 
J0o/ CM^  =  photon  density  of  laser  beam  at  the  scattering  volume 

In  order  to  use  equation  (1),  &  eff  was  determined  by  taking  appropriate  values  of 
Te  and  Ne  and  computing  the  return  to  channels  #1  and  #2  over  the  (sine)2  and 
(cosine) 2  curve.  These  were  summed  and  S  eff  computed  from 

6  eff  =  ^.computer  (2) 

Ne 

This  was  done  for  the  expected  range  of  temperatures  and  concentrations  and 
averaged  to  ±,10%.  The  result  of  these  efforts  produced  a  6  eff.  =  3. 25  x  10"2®cm2 
for  both  the  tunnel  and  arc  experiments.  The  remaining  quantities  in  equation  (1) 
were  then  computed  separately  for  the  arc  and  tunnel  experiments  from  the  geometry 
of  the  experimental  arrangement. 

For  the  arc  runs,  with  0  =  150°,  the  scattering  volume  is  an  ellipitical 
cylinder  bounded  by  the  diameter  of  the  laser  beam  and  the  height  and  width  of  the 
detector  slit.  For  most  of  the  arc  runs  an  entrance  slit  of  50^  x  1. 0  mm  was 
used.  Hence  the  semi-axes  of  the  ellipse  are  1/2  mm  by  2  mm,  the  2  mm  being 
determined  by  the  size  of  the  laser  beam  (  2  mm)  and  0  =  150°.  The  thickness 
of  the  ellipse  being  the  slit  width.  With  these  values  the  scattering  volume  becomes 
V  scat.  =  1.  56  x  10"4  cm2. 
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The  photon  density  was  obtained  by  computing  the  area  of  the  laser  beam 
and  using  the  fact  that  there  are  3.  47  x  1018  0/ joule  at  6943  A.  The  solid  angle 
for  the  system  was  1.23  x  10"^  with  the  resultant  expression, 

Ne  =  0  (measured)  3.  60  x  10*0  (cm-®) 

J  (measured)  ^  ' 

For  the  tunnel  runs,  the  procedure  was  the  same  except  the  working  angle 
0  =  173°  with  the  resultant  expression, 

Ne  -  0  (measured)  6. 80  x  10®  (cm- 3)  (4) 

J  (measured) 

The  J  in  these  equations  must  be  expressed  in  joules. 

The  uncertainty  in  these  working  equations  is  in  the  diverging  of  the  laser 
energy  over  the  laser  beam  area.  In  both  cases  it  was  assumed  that  the  measured 
energy  was  over  a  1.  5  mm  diameter  rather  than  2  mm.  This  was  justified  on 
the  basis  of  a  modified  Gaussian  distribution  of  energy  of  the  laser  beam. 
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APPENDIX  n 


I 

ENGINEERING  DESIGN  DATA  FOR  DETECTOR 
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APPENDIX  m 

COMPUTER  PROGRAM  LISTING  AND  RESULTS 


PROGRAM  NOTES 


Four  subroutines  are  called  which  are  not  herein  listed: 

a.  *  EPMDE  (KODE) 

and  b.  **  CLCINT  (IND,  DX,  EX,  TEMP) 

c.  ***  TNT1  (X,  NTAB,  XTAB,  YTAB,  N,  I  ERR) 

d.  EXIT 

The  purpose  of  (a.)  is  to  suppress  printing  of  error  messages  of  the  form: 

DIV  CHECK  AT  XXXXX  (divide  by  zero) 

EXP  OVERFLO  AT  XXXXX  (number  too  big) 

EXP  UNDERFLO  AT  XXXXX  (number  too  small) 

The  purpose  of  (b. )  is  to  perform  integrations.  Subroutine  (FORTRAN 
function  subroutine)  uses  Simpson's  rule. 

The  purpose  of  (c. )  is  to  interpolate  linearly  by  Lagrange's  formula. 

The  purpose  of  (d. )  is  to  terminate  execution  of  the  program. 


♦Reference:  GE  Space  Div.  'Programmer's  Reference  Manual" 
RM1001,  page  VI-16  • 

♦♦Reference:  GE  Space  Div.  'Engineering  &  Scientific  Computer 
Library"  manual 

♦♦♦Reference:  GE  Information  Systems  "GE-600  Line  FORTRAN  IV", 

CPB-1006G,  Page  70 
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APPENDIX  IV 
ERROR  STATEMENT 


Although  it  is  difficult  to  evaluate  these  test  results  unless  compared  to 
results  obtained  from  a  different  technique,  it  is  instructive  to  compute  the  error 
obtained  in  reading  the  data  and  measurement  of  the  laser  energy.  We  evaluate 
this  for  one  measurement  (Test  #304  of  the  tunnel  experiments)  to  obtain  an  in¬ 
dication  of  the  error.  The  initial  computation  is  the  summation  error  resulting 
in  a  photon  error  ±250  0  on  PM  #1  and  ±590  0  on  PM  #2.  It  is  noted  that  these 
values  are  greater  than  the  statistical  error.  Hence  PM  #1  is  27, 500  ±250  and 
PM  #2  is  16, 500  ±590.  The  summation  of  these  two  channels  results  on  an  error, 

£l)2  =  44, 000  ±  \/  (250)2  +  (590)2  (1) 

or 

2^  44» 000  ±  650  0  (2) 

To  obtain  Ne,  we  divide  this  value  by  the  energy  of  the  laser  beam  assumed  to  be 
accurate  to  ±10%  or  0. 11  ±  .  01.  The  probable  error  can  be  obtained  from, 


P.E.r  }or  BA 

where  a  is  error  in  A  and  b  is  the  error  in  B.  Substituting  we  obtain. 


p-  E.  =  *  -L (4.4*10^  IQ"1  (*) 
0A“  (WXIo'1)* 

Evaluating  we  obtain, 

P.  E.  =  ±  3.  6  x  10^  0/joule  (4) 
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Since  the  total  signal  is  (4.  4  x  10^)/0. 11  or  4.  0  x  10®  0/ joule  and  Ne  is 
proportional  to  this  value  we  conclude  that  the  summation  process  yields 
a  measurement  of  ±10%. 

To  obtain  the  electron  temperature,  we  ratio  PM  #2  output  to  PM  #1 
output  or 


jO  =  16,  500  ±  590  (5) 

/  27,  500  ±  2F0 

Evaluating  this  expression  in  the  same  way,  we  obtain 
yO  =  0.  600  ±  0. 02 

Referring  to  FIG  19,  at  6000°K  and  with  Nesi  2. 8  x  10^6  this  corresponds  to 
an  error  in  temperature  of  ±200°K.  A  more  accurate  measurement  of  tem¬ 
perature  can  probably  be  obtained  by  more  careful  measurement  of  the  voltage. 
However,  the  measurement  of  electron  concentration  is  set  by  the  accuracy  of 
the  laser  energy  monitor. 


t 


74 


AEDC-TR-71-69 


APPENDIX  V 


REPRODUCTION  OF  PAPER  BY  C.  J  HARRIS 


AEDC-TR-71-69 


CORRELATION  OF  INVISCID  AIR  NONEQUILIBRIUM  SHOCK  LAYER  PROPERTIES* 


Clarence  J.  Harris,  Staff  Scientist 
General  Ele'ctrlc  Company 
Valley  Forge,  Pa. 


Abstract 


It  has  already  been  shown  that  the  entropy  value  de¬ 
fined  under  equilibrium  stagnation  conditions  in  the 
plenum  or  reservoir  region  of  a  nozzle  readily  corre¬ 
lates  and  permits  the  prediction' of  the  quantitative  non¬ 
equilibrium  chemical  and  electron  composition  of  the 
downstream  converging-diverging  nozzle  expanding  air¬ 
flow.  This  correlation  based  on  plenum  entropy  also  can 
be  applied  to  define  the  axial  station  in  a  given  nozzle 
where  transition  from  equilibrium  to  "frozen"  non- 
equilibrium  flow  occurs.  In  this  paper  an  analog  is  de¬ 
veloped  relating  this  "entropy  correlation"  for  the 
axisymmetrlc  nozzle  situation  to  the  shock  layer  flow 
situation  about  a  conical  body  with  nose  bluntness.  The 
analog  as  initially  developed  was  limited  to  considering 
the  innermost  stream  tube  of  an  lnviscld  nonequilibrium 
shock  layer  airflow.  The  earlier  nozzle  entropy  corre¬ 
lation  was  used  to  predict  the  frozen  nonequilibrium 
properties  in  this  stream  tube  and  these  predicted  prop¬ 
erties  then  compared  to  shock  layer  values  derived  from 
a  detailed  nonequilibrium  flow  field  numerical  calcula¬ 
tion.  The  resulting  comparison  indicates  that  the 
"nozzle-shock  layer  analog"  and  entropy  correlation  are 
useful  in  predicting  shock  layer  properties  thus  reducing 
the  need  for  extensive  numerical  calculations.  Additional 
applications  and  limitations  of  the  newly  developed 
"analog"  technique  are  evaluated,  particularly  as  they 
aid  in  quickly  defining,  as  a  function  of  vehicle  altitude 
and  flight  velocity,  the  properties  at  the  inviscld  shock 
layer  flow>-vlscous  boundary  layer  edge.  Preliminary 
experimental  laser  plasma  probe  measurements  also  are 
presented  of  electron  density/temperature  in  a  high 
enthalpy  shock  tunnel  nozzle  air  flow.  These  results 
suggest  an  entropy  correlation  for  predicting  non- 
equilibrium  electron  temperatures  in  frozen  heavy  parti¬ 
cle  flows. 


Nomenclature 


(b) 

». 


h 

H/RT0 

i 

Nes 

P5 

rl 

rb 

r„ 


Annular  stream  tube  area  at  Station  b 
(see  Fig.  4) 

Annular  stream  tube  area  at  sonic  point 
Area  ratio 

Bluntness  ratio,  r_/r, 
n  b 

Altitude 

Enthalpy  ratio 

Nozzle  scale  parameter 

Free  stream  Mach  number 

Electron  concentration  at  stagnation  point 

Plenum  or  reservoir  stagnation  pressure 

Centerbody  radius  at  sonic  point 

Centerbody  radius  at  Station  b 

Nose  radius 


S/R 

t5 


Entropy  ratio 

Plenum  or  reservoir  stagnation  temperature 
Free  stream  velocity 

Thickness  of  Innermost  inviscid  stream  tube 


I.  Introduction 


The  phenomena  associated  with  the  expansion  of  air 
in  a  nozzle  to  hypersonic  velocities  have  been  studied  by 
several  authors.  Early  analytical  work  on  flow  expan¬ 
sions  and  nonequilibrium  effects  was  carried  out  by 
Brayt*-)®  on  monoatomlc  gases.  More  recently  Lord! 
and  Mates  (3)  did  a  detailed  computer  program  study  for 
high  enthalpy  air  flow*  v* tv  nonequilibrium  effects 
present.  Harris^)  utilized  the  results  of  Reference  3 
and  similar  computer  calculations  by  Boyer^5)  and 
Eschenroeder  et  al^  to  show  that  the  nonequilibrium 
chemical  composition  in  the  expanding  nozzle  airflow 
could  be  readily  correlated  as  a  function  of  plenum  or 
reservoir  entropy.  Harris  et  al(7)  also  showed  that  a 
similar  "entropy  correlation”  was  valid  for  specifying 
the  nonequilibrium  electron  concentration  values  in  an 
expanding  air  flow.  Both  Harney^8)  and  Ring  and 
Johnson^  have  made  extensions  to  the  entropy  corre¬ 
lation  concept  offered  in  References  4  and  10.  This 
simplified  entropy  correlation  for  handling  energetic 
nozzle  air  flows  has  had  direct  application  in  under¬ 
standing  and  defining  hypersonic  and  hypervelocity  test 

facilities  performancef10H1l)(12)(13)<  in  Ujis  paper  the 
nozzle  nonequilibrium  entropy  correlation  approach  is 
adapted  to  the  case  of  hypersonic  shock  layer  flow 
around  a  blunt,  slender  axisymmetrlc  body  and  accord¬ 
ingly  used  to  predict  certain  pertinent  properties  in  the 
shock  layer.  This  work  is  discussed  in  greater  detail 
in  Reference  20. 


II.  Analysis 


'■'Nozzle-Shock  Layer  Analog" 


A.  Nozzle 


The  general  conditions  required  for  the  valid  appli¬ 
cation  of  the  entropy  correlation^10)  are: 

1.  The  plenum  or  stagnated  reservoir  reacting 
gas  (air)  is  in  equilibrium. 

2.  The  reacting  flow  stays  in  equilibrium  in 
expanding  from  the  plenum  to  the  nozzle 
throat  (sonic  point). 

3.  Viscous  effects  are  ignored. 
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This  is  illustrated  using  the  air  Mollier  Diagram  in  Fig¬ 
ure  1.  The  "departure  from  equilibrium"  boundary 
(Figure  1)  was  defined  in  Reference  3  and  the  "frozen 
flow"  limit  on  the  Air  Mollier  Diagram  is  as  defined  in 
References  10  and  11.  The  entropy  correlation  has 
meaningful  application  only  when  the  plenum  operating 
points  (such  as  A,  B,C,D,E  and  F  in  Figure  1)  are  above 
the  departure  from  equilibrium  boundary  on  the  Mollier 
Diagram.  Should  the  plenum  operating  point  fall  below 
this  boundary  or  the  frozen  flow  limit,  the  entropy  cor¬ 
relation  is  not  applicable. 

The  Air  Mollier  Diagram  accurately  defines  equi¬ 
librium  thermodynamic  and  kinetic  properties  (tem¬ 
perature,  pressure,  density,  velocity,  and  Z)  above  the 
departure  from  equilibrium  boundary.  Below  this 
boundary  and  the  frozen  flow  line,  finite  rate  chemistry 
occurs  and  the  equilibrium  thermodynamic  and  kinetic 
properties  are  no  longer  valid.  However,  it  has  been 
shown^W®)^®)  that  for  an  expanding  flow  the  equilibrium 
and  nonequilibrium  particle  velocity  and  density  differ 
very  tittle. 

With  this  kind  of  a  Mollier  Diagram  representation 
of  the  nozzle  flow  process  it  is  also  possible  to  define 
the  area  ratio  (A/A*)  value  at  which  "freezing"  occurs. 
The  performance  ljmit  line  for  the  GE  High  Enthalpy 
Shock  Tunnel'1  'is  also  included  in  Figure  1  and 
provides  performance  equal  to  or  in  excess  of  similar 
facilities^®).  So,  it  may  be  expected  that  in  most  cur¬ 
rently  operating  test  facilities  nonequilibrium  freezing 
occurs  at  small  area  ratio  values  or  just  downstream  of 
the  nozzle  throat.  The  original  entropy  correlation  of 
frozen  nonequilibrium  chemical  (Figure  2)  and  electron 
(Figure  3)*  composition  data  (References  4,  10,  and  7) 
was  carried  out  at  and  presented  for  A/ A*  =  10^.  None¬ 
theless,  these  correlations  are  still  valid  for  the  much 
smaller  A/A*  values  such  as  defined  under  the  condi¬ 
tions  in  Figure  1.  This  area  ratio  concept  for  defining 
the  station  where  nonequilibrium  frozen  air  flow  occurs 
is  applied  to  the  shock  layer  flow  situation. 

B.  Shock  Layer 

This  section  deals  with  establishing  an  analog  be¬ 
tween  nozzle  flow  and  shock  layer  flow  and,  accordingly, 
applies  to  the  shock  layer  flow  situation  and  nozzle  pre¬ 
diction  techniques  discussed.  A  rigorous  mathematical 
solution  to  the  equilibrium  shock  layer  flow  field  prob¬ 
lem  was  derived  quite  early  by  Gravalos  and  his  fellow 
workers^1®).  Gravalosf1^)  has  recently  summarized  the 
analytical  work  in  this  field  for  the  flow  of  inviscid,  re¬ 
acting  gas  mixtures  about  axisymmctric  bodies. 

The  typical  shock  layer  flow  situation  to  be  con¬ 
sidered  is  given  in  Figure  4.  The  center  body  is  a  cone 
with  a  hemispherical  nose.  Thus  the  stream  tube 
crossing  the  bow  shock  wave  at  (a)  is  eventually  stag¬ 
nated  at  (s)  and  then  expands  to  sonic  velocity  at  (*)  and 
becomes  supersonic  downstream  of  this  point.  Thus, 
region  (s)  in  the  shock  layer  flow  is  analogous  to  the 
plenum  region  in  the  nozzle  flow  situation  (see  Figure  4a) 


and  (*)  in  the  shock  layer  situation  is  analogous  to  the 
throat  in  the  nozzle  flow  case  (see  Figure  4b. ).  Further, 
the  supersonic  expansion  of  the  center  body  stream  tube 
downstream  of  the  sonic  point  is  analogous  to  the  super¬ 
sonic  expansion  in  the  nozzle  case.  Therefore,  the 
present  analysis  is  restricted  to  an  evaluation  of  the 
center  body  stream  tube  and  regions  such  as  (s),  (*), 
and  later  on  (b),  in  the  shock  layer  flow.  It  is  necessary 
to  require  that  the  reacting  stagnated  air  in  region  (s)be 
in  equilibrium  and  that  it  undergo  an  equilibrium  expan¬ 
sion  to  the  sonic  point  (*).  It  is  next  desirable  to  define 
an  equivalent  area  ratio  A(b)/A(*)  for  the  center  body 
stream  tube  at  and  beyond  which  the  chemistry  should  be 
"frozen"  in  its  nonequilibrium  state.  This  equivalent 
area  ratio  is  defined  using  the  following  assumptions: 

1.  The  total  shock  layer  is  inviscid 

2.  The  whole  reacting  airflow  region  between  the 
bow  shock  and  body  and  from  the  stagnation 
line  to  the  sonic  line  is  in  equilibrium  and 

3.  The  center  body  stream  tube  thickness,  x«r 
(Figure  4). 

Under  these  first  two  conditions,  Edsall'  showed  for 
a  wide  range  of  free  stream  Mach  numbers  (M  =  3.  5  to 
35)  that 

r  *■  (0.6  to  0.7)  rn  (1) 


Therefore,  the  annular  stream  tube  area  at  (*)  and  (b) 
will  be 


A(*)  =  (X  +  rj)2  -  n  rx2 

(2) 

A(b)  =  ff(XTrb)2-ffrb2 

(3) 

then 

.  2 

A.  /A  „  =  x  +  2xr.  ,  2  „  . 

(b)  (*)  b/(x  +2x  r^) 

(4) 

since 

B  (bluntness  ratio)  =  r  , 
n/rb 

(5) 

and  r„  *»  0.  65  r 

1  n 

(6) 

then 

x  +  2r 

.  /.  n/B 

A(b)/A(*)fts  x+1.3rn 

(7) 

let 

i  =  x/r„  =  f  (M».  h' 

(8) 

*Note  that  these  results  supercede  those  offered  in  Reference  7  where  the  data  had  been  mistakenly  miss-plotted. 
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then 

A  /A  M  f  +  2/B 
(b)  (*)  C  +  1.3 

for  §  «  1  and  B  <  1 

W2/UB 


(9) 


(10) 


Assumption  No.  2  was  necessary  for  arriving  at  a 
general  definition  as  to  the  location  of  the  sonic  point  on 
the  center  body  (i.e. ,  Eq.  6).  The  next  step  is  to  deter¬ 
mine  the  stagnation  enthalpy  and  entropy  values  at  (s)  for 
the  flight  conditions  under  consideration.  If  this  opera¬ 
ting  point  falls  above  the  departure  from  equilibrium 
boundary  on  the  Air  Mollier  Diagram  (Figure  1),  then  the 
equivalent  area  ratio  or  bluntneBs  ratio  (B)  at  which 
freezing  will  occur  may  be  determined.  Once  this  center 
body  station  (or  area  ratio)  is  defined,  the  chemical  com¬ 
position  and  electron  concentration  in  the  center  body 
stream  tube  at  points  downstream  of  this  station  may  be 
predicted  (from  Figure  2  and  Figure  3,  respectively) 
based  on  the  stagnation  point  entropy  value. 


composition  for  the  center  body  stream  tube  at  this  body 
Btation  is  compared  for  the  four  flight  conditions  (in  Fig¬ 
ure  5)  with  the  composition  based  on  the  nozzle  entropy 
correlation^10).  The  species  considered  are  02»  N2. 
O,  and  N.  The  symbols  are  for  the  numerically  calcu¬ 
lated  Bhock  layer  values  and  the  solid  lines  are  from  the 
nozzle  correlation^4^10). 

With  the  exception  of  the  order  of  magnitude  dis¬ 
agreement  between  the  calculated  value  of  N  and  the  value 
for  N  predicted  by  the  correlation  at  the  50K  feet  condi-  _ 
tion,  the  overall  agreement  is  quite  good  and  the  corre¬ 
lation  appears  to  predict  the  chemical  properties  in  the 
frozen  nonequilibrium  portion  of  the  inviscid  center  body 
stream  tube  quite  well.  A  similar  comparison  is  offered 
in  Figure  6  for  the  electron  concentration.  While  the 
general  quantitative  level  of  the  numerically  calculated 
shock  layer  results  follows  that  of  nozzle  correlation,  the 
trend  and  agreement  between  the  two  is  not  clear.  Since 
the  prediction  of  shock  layer  electron  concentration  as  a 
function  of  flight  trajectory  is  of  importance  to  many 
investigators,  the  results  of  Figure  6  are  evaluated  in 
more  detail  in  the  next  section. 


Comparison  with  Numerical  Calculations 


111.  Applications  and  Limitations 


Inviscid,  nonequilibrium  sphere-cone  body  shock 
layer  calculations  have  recently  been  performed  by 
■McMenamin  and  O'Brien^19).  The  results  from  these 
calculations  were  compared  to  the  results  based  on  the 
prediction  techniques  discussed  above.  While  the  results 
of  the  calculations^19)  are  quite  extensive,  covering  the 
full  shock  layer  flow  field  and  also  considering  the  prob¬ 
lem  of  inviscid  flow-viscous  boundary  layer  interaction, 
attention  here  is  again  restricted  only  to  the  center  body 
stream  tube  for  the  case  of  purely  inviscid  flow;  i.e.  , 
where  no  boundary  layer  is  present  (although  the  center- 
body  shape  may  have  been  altered  to  account  for  the 
boundary  layer  displacement  thickness,  6*). 

The  center  body  configuration  selected  was  one  which 
had  a  cone  angle  of  <10°  and  slight  cone  bluntness.  In¬ 
viscid  flow  field  results  were  chosen  for  four  flight  con¬ 
ditions.  The  prediction  techniques  proposed  herein  were 
compared  to  the  results  from  these  flight  conditions.  The 
four  flight  conditions  are: 


h  (ft) 

wm  (ft/ sec) 

1. 

50  x  103 

19,200 

2. 

100  x  103 

21,500 

3. 

150  x  103 

21,950 

4. 

200  x  103 

22,000 

Detailed  flow  field  calculations^9)  were  available  for 
three  normals  along  the  center  body  surface  for  each  of 
these  flight  conditions.  The  most  downstream  of  these 
normals  is  at  a  body  station  yielding  an  equivalent  area 

ratio  value  of  approximately  50,  thus  well  into  the 
"frozen”  nonequilibrium  flow  regime  (20)  on  the  body. 
The  numerically  calculated^19)  shock  layer  chemical 


Flight  Trajectory  and  Vehicle  Size 

The' foregoing  analysis  has  been  restricted  to  the 
case  of  a  fully  inviscid  shock  layer  flow  and  the  predic¬ 
tion  techniques  have  been  compared  only  to  one  particular 
body  shape  and  size.  The  first  restriction  requires  that 
the  prediction  techniques  should  only  be  applicable  in  the 
flight  regime  where  the  inviscid  shock  layer  flow  situ¬ 
ation  is  approximated.  Probstein(21)  has  defined  seven 
regimes  that  develop  during  atmospheric  re-entry  for  the 
flow  field  environment  about  a  body  having  a  spherical 
stagnation  regime.  For  a  nose  radius  of  0. 5  inch,  the 
regime  of  interest,  the  viscous  boundary  layer-inviscid 
shock  layer  regime,  occurs  at  altitudes  of  175K  feet  and 
below.  For  a  nose  radius  of  1  foot  the  maximum  altitude 
would  be  adjusted  to  240K  feet. 

The  variation  in  nozzle  size  or  "time -rate -of-  ■ 
expansion”  and  its  effect  on  the  entropy  correlation  is 
considered  in  detail  in  References  8  and  9.  In  the  corre-. 
lation  of  References  4,  7,  and  10  limited  consideration 
was  given  to  the  physical  scale  of  the  nozzle  and  its  effect 
on  the  entropy  correlation.  Figure  2  shows  that  a  scale 
factor  (l )  variation  of  4.  7  does  not  greatly  affect  the 
nozzle  entropy  correlation  as  applied  to  chemical  compo¬ 
sition.  Ring  and  Johnson^9)  looked  at  a  wide  range  of 
values  of  the  nozzle  size  and  shape  similarity  parameter 
(Tc)  and  their  influence  on  the  entropy  correlation.  Fig¬ 
ures  3  and  6  indicate  that  the  dimensional  scale  factor 
does  have  a  pronounced  effect  on  the  nozzle  entropy  cor¬ 
relation  as  regards  electron  concentrations. 

Further,  Eschenroeder(22)  contends  that  ionization 
freezing  in  a  nozzle  takes  place  upstream  of  neutral 
chemical  freezing.  So  while  the  entropy  correlation  and 
nozzle -shock  layer  analog  may  have  direct  application  in 
the  case  of  neutral  chemical  composition,  more  detailed 
consideration  must  be  given  to  the  electron  concentration 
situation.  This  will  be  discussed  further  when  the 
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experimental  laser  probe  measurements  are  presented 
of  electron  density/temperature.  Reference  22  also 
points  up  the  differences  in  electron  concentration  values 
that  result  from  considering  a  "simple  kinetic"  model  as 
compared  to  a  "complex  kinetic"  model.  The  simple 
kinetic  model  allows  for  only  nitric  oxide  ionization. 

The  complex  kinetic  model  allows  for  the  additional 
species  of  0+,  N+,  and  N2^.  For  this  reason  the  com¬ 
parison  of  the  nozzle  entropy  correlation  of  electron 
concentration  to  the  shock  layer  numerical  calculations 
(Figure  6)  is  replotted  in  Figure  7  to  account  for  this. 

Reference  20  showed  that  for  both  the  nozzle  and 
shock  layer  cases  the  ionized  species  are  determined 
using  the  nitric  oxide  dissociative-recombination  re¬ 
action  only.  Eschenroeder(22)  argues  that  the  model  in¬ 
volving  this  reaction  and  the  three  body  NO  reaction  is 
correct  for  nozzle  flows  where  plenum  pressures  are 
usually  quite  high.  But,  for  shock  layer  flow  where  the 
stagnation  pressure  is  much  lower,  die  complex  kinetic 
model  is  applicable.  The  pressure  argument  may  in  fact 
really  be  an  entropy  argument.  That  is,  below  a  certain 
plenum  entropy  value  the  simple  kinetic  model  is  valid 
and  above  that  entropy  value  the  complex  kinetic  model 
holds.  This  can  be  argued  from  Figures  2  or  5.  That 
is,  for  entropy  values  up  to  S/R  =  38,  the  NO  concentra¬ 
tion  exceeds  the  N  atom  concentration  and  electron  pro¬ 
duction  might  be  expected  to  be  controlled  mainly  by  the 
three  body,  M  +  NO+  +  e”5  NO  +  M  reactions  and  par¬ 
tially  by  the  dissociative-recombination  reaction 
NO+  +  e"SN  +  0<28>.  Above  S/R  =  38  the  latter  reaction 
should  be  a  major  electron  production  contributor  and  the 
N+  reaction  becoming  quite  an  effective  contributor, 
along  with  the  0+  and  Ng+  reactions. 

Thus ,  the  complex  kinetic  model  appears  to  be  valid 
at  the  higher  entropy  values  and  the  simple  kinetic  model 
valid  at  the  lower  S/R  values.  This  is  a  postulation  and 
requires  further  investigation.  Nevertheless ,  the  infor¬ 
mation  in  Figure  7  was  replotted  without  modifying  the 
nozzle  results  and  the  two  shock  layer  values  at  the  lower 
entropy  conditions  (50K  feet  and  100K  feet  altitude).  The 
150K  feet  altitude  shock  layer  value  was  replaced  by  its 
equilibrium  value  since  it  is  closer  to  the  complex  kinetic 
model  nonequilibrium  value *22*  and  the  nonequilibrium 
complex  kinetic  model  value  obtained  in  Reference  22 
was  included. 

Using  these  values  for  the  shock  layer  flow  where 
the  stagnation  point  entropy  value  is  large  causes  the 
data  to  lie  an  a  relatively  smooth  curve.  Also  presented 
in  Figure  7  are  experimental  data  by  Dunn  and 
Lordi^23, 28)  obtained  in  recent  shock  tunnel  nozzle 
studies.  The  data  are  from  an  experiment  in  which 
l-z  3.26  cm.  These  experimental  results  indicate  that 
the  chemical  kinetic-thermodynamic  model  of  Refer¬ 
ences  3,5,  and  6  slightly  overpredicts  the  measured 
frozen  nonequilibrium  electron  concentration  in  the 
nozzle.  These  data  and  electron  density/temperature 
data  obtained  in  this  study  are  discussed  later  in  this 
paper. 

The  present  nozzle-shock  layer  analog  specifies  no 
particular  restrictions  on  center  body  cone  half  angle. 
Figure  1  and  the  Equations  (1)  to  (10)  imply  only  a  blunt¬ 
ness  ratio  dependence.  Therefore,  the  analog  should 


hold  for  a  range  of  center  body  cone  half  angles .  How¬ 
ever,  the  restrictions  on  this  argument  need  to  be  more 
thoroughly  developed  and  investigated: 

Application  to  Flight 

An  additional  specification  discussed  earlier  was  the 
requirement  that  the  stagnation  enthalpy-entropy  value  of 
the  air  be  such  as  to  place  the  equivalent  plenum  opera¬ 
ting  point  ((s)  Figure  4)  above  the  "departure  from  equi¬ 
librium"  boundary  and  the  "frozen  flow"  line  on  the  Air 
Mollier  Diagram  (Figure  1).  The  frozen  flow  line  and 
the  maximum  flight  altitude  at  which  an  inviscid  shock 
layer-viscous  boundary  layer  flow  situation^2*)  is  ap¬ 
proximated  may  be  used  jointly  to  define  the  limits  of 
applicability  of  the  entropy  correlation  in  predicting 
shock  layer  properties.  These  limits  are  shown  in  Fig¬ 
ures  8  and  9.  Figure  8  presents  the  entropy  correlation 
of  shock  layer  flow  properties  as  a  function  of  vehicle 
flight  conditions.  Flight  velocity  (vro)  and  center  body 
stream  tube  stagnation  point  enthalpy  (HS/RT0)  are  pre¬ 
sented  versus  entropy  (S/R).  Lines  of  constant  flight 
altitude  (h)  and  stagnation  point  electron  density 
(Ne8(equiiib.  are  a^so  deluded.  The  four  flight  condi¬ 
tions  that  were  used  are  shown  (filled  circles)  as  well  as 
the  flight  point  from  Reference  22  (open  circle).  The 
altitude  boundaries  for  achieving  a  viscous  boundary 
layer-inviscid  shock  layer  flow  are  presented  for  both  a 
0.  5  inch  and  1.0  foot  nose  radius  center  body. 

The  other  limit  line  is  the  frozen  flow  line.  In  Ref¬ 
erences  4  and  10  this  line  was  defined  based  upon  the 
degree  of  oxygen  dissociation  (in  upper  Figure  8  denoted 
by  hatching  directed  towards  the  left  on  the  frozen  flow 
limit  line).  The  percentage  of  undissociated  oxygen 
derivable  from  the  entropy  correlation  is  tabulated  along 
the  entropy  (S/R)  axis  of  upper  Figure  8.  The  frozen 
flow  limit  line  was  extended  by  considering  the  total  de¬ 
gree  of  air  dissociation  and  accounting  for  the  nitrogen 
dissociation  (in  upper  Figure  8  denoted  by  hatching 
directed  towards  the  right  on  the  frozen  flow  limit  line). 
Therefore,  in  the  region  in  Figure  8  that  is  below  175K 
feet  for  a  0. 5  inch  nose  radius  and  which  is  above  the 
frozen  flow  line,  it  is  possible  to  apply  the  entropy  cor¬ 
relation  in  predicting  the  center  body  stream-tube  shock 
layer  flow  properties.  This  can  be  done  in  the  acceptablf 
regime  at  any  point  along  the  flight  trajectory  of  a 
slightly  blunted  vehicle  since  the  stagnation  point  entropy 
is  specified  by  altitude  and  velocity.  Other  limits  are 
those  in  which  the  chemical  model  is  valid.  The 
model^H24)  is  taken  to  be  valid  below  200K  feet  and  for 
less  than  30(17). 

Presented  below  the  flight  condition  plot  in  Figure  8 
is  a  replot  of  the  entropy  correlation  of  neutral  species 
and  electron  concentration.  One  need  only  to  use  the  top 
half  of  Figure  8  to  determine  the  stagnation  point  entropy 
value  at  any  given ‘point  on  a  flight  trajectory  (Vg,  and  h). 
Then,  by  projecting  straight  down  along  this  entropy 
value  to  the  curves  below,  an  estimate  of  the  electron 
concentration  is  obtained  as  well  as  fairly  accurate 
quantitative  values  of  the  neutral  chemical  species  con¬ 
centration.  Thus,  for  example  a  vehicle  at  100K  feet 
altitude  and  flying  at  25, 000  ft/sec  is  predicted  to  have 
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a  downstream  center  body  stream  tube  electron  density 
value  of  3  x  10  (mol/gm  mix. )  and  the  neutral  chemical 
species  composition  specified  at  the  entropy  (S/R)  value 
of  46.  The  stagnation  point  equilibrium  electron  density 
value  will  be  10*®  electrons/cm®. 

Boundary  Layer  Effects 

The  assumption  made  so  far  of  a  fully  inviscid  shock 
layer  flow  is  only  an  approximation  of  the  real  shock 
layer  flow  situation  which  Involves  not  only  an  inviscid 
flow  regime  but  also  a  viscous  boundary  layer.  It  was 
the  need  to  treat  this  realistic  flow  situation  that  led  to 
the  analytical  study  and  numerical  calculations  of  Ref¬ 
erence  19.  In  that  study  detailed  consideration  was 
given  to  the  mixing  or  "swallowing"  that  takes  place 
between  the  Inviscid  flow  and  the  viscous  boundary  layer 
and  the  resulting  effect  on  the  properties  in  the  full 
shock  layer.  In  this  investigation  consideration  was 
given  to  the  influence  on  electron  concentration  of 
catalytic  wall  -  noncatalytic  wall  effects ,  laminar 
boundaiy  layer,  and  impurities  in  the  boundary  layer. 
This  discussion  is  given  in  Reference  20, 

In  general,. the  inviscid  solutions  yield  Ng  values  an 
order  of  magnitude  greater  than  those  obtained  with  a 
pure  air  laminar  boundary  layer  and  a  catalytic  wall. 
However,  the  measurements  discussed  in  the  next  sec¬ 
tion  indicate  that  adjustments  need  to  be  made  in  the 
inviscid  theoretical  model  that  accordingly  lead  to  a 
factor  of  three  or  four  reduction  in  the  predicted  electron 
concentration. 

Since  the  analog  is  based  on  the  entropy  correlation 
that  in  turn  assumes  negligible  changes  in  entropy 
during  the  expansion  process  <29),  the  shock  layer  pre¬ 
dictions  are  only  goocf  along  the  body  from  the  upstream 
"freeze"  station  to  the  downstream  station  where  pro¬ 
nounced  entropy  changes  begin  to  occur.  Downstream  of 
this  station  the  prediction  techniques  are  only  good  for 
making  estimates. 

IV.  Experimental  Study 

Laser  Plasma  Probe-Shock  Tunnel  Electron  Density/ 
Temperature  Measurements 

Limited  data  were  obtained  of  the  electron  density 
and  temperature  in  the  expanding  air  flow  of  the  GE  High 
Enthalpy  Shock  Tunnel  A  recently  developed 

Laser  Plasma  Probe^27^  was  used  to  obtain  these  data 
at  a  A/A*  value  of  3.6  where  the  flow  is  predicted  to  be 
in  equilibrium  and  of  high  density  <p/p0  ~  1).  The 
measured  temperature  values  agreed  within  10%  with 
the  theoretical  equilibrium  heavy  particle  gas  temper¬ 
ature  values.  The  measured  electron  density  values 
were  within  20%  of  the  theoretical  equilibrium  values. 
These  preliminary  results  Imply  that  the  electrons  and 
heavy  particles  are  in  "equilibrium"  with  each  other  In 
the  theoretically  predicted  equilibrium  expanded  flow 
regime.  From  the  entropy  correlation,  these  experi¬ 
mental  results,  and  the  only  other  experimental  work 
of  this  type  available,  that  of  Dunn  and  Lordl  (23)(28)t 
a  model  for  the  "nonequilibrium"  electron  temperature 


characteristics  was  postulated.  It  is  postulated  that  the 
electron  temperature  and  heavy  particle  temperature 
are  the  same  in  the  equilibrium  expanding  flow  regime. 
When  the  flow  expands  to  its  "freeze"  point  the  electron 
and  heavy  particle  temperatures  uncouple.  For  an 
example  look  at  the  plenum  or  reservoir  point  E  in 
Figure  1.  The  flow  expands  in  equilibrium  from  this 
point  to  an  area  ratio,  A/A* ,  of  4.  At  this  expansion 
value,  frozen, nonequilibrium  flow  exists.  Further,  the 
decay  in  electron  temperature  with  increasing  expan¬ 
sion  lags  the  heavy  particle  temperature  decay.  This 
is  illustrated  in  Figure  10,  using  the  results  of  Refer¬ 
ences  23,  27,  28.  The  postulated  "electron  initial  lag 
temperature"  is  uniquely  defined  by  the  entropy 
correlation  and  would  be  as  presented  in  Figure  11. 
Verification  of  this  postulation  may  be  accomplished  by 
an  experimental  study  in  which  the  electron  density  and 
temperatures  are  measured  both  in  the  equilibrium  and 
nonequilibrium  expanded  flow  regimes.  At  present  the 
laser  plasma  probe  appears  to  be  the  better  technique 
for  accomplishing  this.  Other  techniques  such  as 
Langmuir  probes,  time-of-Qight  mass  spectrography, 
and  microwave  probes  are  severely  limited  as  to  the 
nozzle  flow  regimes  In  which  they  are  applicable. 

V.  Conclusions 

Within  the  appropriate  flight  regime  as  defined  in 
Figures  8  and  9,  the  entropy  correlation  and  nozzle- 
shock  layer  analog  may  be  used  to  predict  the  frozen 
nonequillbrium  properties  in  the  inviscid  inner  shock 
layer  stream  tube.  Comparison  to  an  analysis  where 
the  viscous  boundary  layer  Is  accounted  for  Bhows  that 
this  inviscid  assumption  when  corrected  on  the  basis  of 
recent  experimental  data  closely  predicts  the  electron 
concentration  value  in  the  body  stream  tube<20\ 

Mapy  facets  of  the  proposed  shock  layer  prediction 
technique  need  to  be  more  critically  evaluated.  These 
include  the  effect  of  center  body  physical  scale,  the 
effect  of  center  body  cone  angle,  and  the  correctness  of 
the  thermodynamic-kinetic  air  model  presently  being 
used.  Dunn  and  Lordl 's  experimental  results'22^  and 
the  results  of  the  preliminary  laser  plasma  probe 
Study  reported  herein  Imply  a  need  to  modify  the 
ionized  species  model.  The  present  prediction  study 
also  implies  that  the  influence  of  the  various  re¬ 
actions  suggested  in  the  ionized  species  model  may  be 
more  directly  dependent  upon  entropy  than  upon  pressure 
as  previously  suggested  by  Eschenroeder^22'  and  that 
the  nonequilibrium  electron  temperature  character¬ 
istics  may  be  specified  by  an  entropy  correlation. 

The  entropy  correlation  and  nozzle-shock  layer  analog 
prediction  technique  offers  several  worthwhile  advan¬ 
tages.  It  provides  a  great  reduction  in  the  machine 
computation  required  in  carrying  out  point -by-point 
flight  trajectory  studies  to  define  or  estimate  important 
properties  such  as  electron  concentration  and  neutral 
chemical  composition  in  the  shock  layer.  And  it  also 
implies  the  possibility  of  using  nozzle  flows  in  experi¬ 
mental  studies  to  simulate  some  of  the  Important 
phenomena  taking  place  in  certain  important  portions  of 
center  body  shock  layer  Qowt2®).  Also  see  Ref.  30. 
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The  applicability  of  the  shock  layer  prediction 
techniques  offered  is  as  shown  in  Figure  9.  The 
Mg,  =  30  and  200K  ft  altitude  limits  define  the  regime 
below  which  the  chemical  model  is  applicable  (17)  (24). 

In  addition,  the  maximum  altitude  limit  is  further  defined 
by  the  need  for  a  boundary  layer-lnviscid  shock  layer 
flow  situation  and  thus  is  a  function  of  centerbody  nose 
radius,  rn(21\  The  other  limit  is  that  below  which 
equilibrium  reacting  chemistry'  is  no  longer  achieved  at 
the  stagnation  point.  Therefore,  the  prediction  technique 
is  generally  applicable  from  200K  ft  for  v«  between  22 
to  33K  ft/ sec  to  h  *>  50K  ft  and  below  for  between  12K 
ft/sec  to  30K  ft/second. 
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Figure  1.  Air  Molller  Diagram  Showing 
Equilibrium  and  Nonequilibrium  Regions 
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Figure  2.  Nonequilibrium  Chemical  Concentrations  vs 
Entropy  for  Expanding  High  Enthalpy  Air  Flows 
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Figure  3.  Nonequilibrium  Electron  Concentration  Correlation  with  Reservoir  Entropy  at  A/A*  =  10' 
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a.  SHOCK  LAYER  FLOW  SITUATION  (INVISCID) 


b.  COMPARABLE  NOZZLE  SITUATION 
Figure  4.  Nozzle  Flow-Shock  Layer  Analog 
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Figure  9.  Flight  Regime  in  which  Prediction  Technique  la  Applicable 


Figure  10.  Composite  of  Electron  Temperature  Data  in  Shock  Tunnels 
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Figure  11.  Postulated  Conditions  at  which  Electron  Temperature 
and  Heavy  Particle  Temperature  Begin  to  Differ 
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Thomson  scattering  radiation  has  been  measured  in  terms  of 
electron  concentration  and  temperature  in  a  free  jet  facility  and  a 
shock  tunnel  using  a  1/2  joule,  Q-switched  ruby  laser.  A  detector 
system  is  described  which  attenuates  the  primary  ruby  wavelength  by 
10  ,  discriminates  between  orthogonal  polarization  states,  analyzes 
the  radiation  spectrally  and  provides  a  measure  of  radiation  background. 
The  scattering  volume  was  approximately  3  x  10“4  cm^  for  the  jet  fa¬ 
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